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Abstract 

Realization of graphene moire superstructures on the surface of 4d and 5d transition 
metals offers templates with periodically modulated electron density, which is respon- 
sible for a number of fascinating effects, including the formation of quantum dots and 
the site selective adsorption of organic molecules or metal clusters on graphene. Here, 
applying the combination of scanning probe microscopy /spectroscopy and the density 
functional theory calculations, we gain a profound insight into the electronic and 
topographic contributions to the imaging contrast of the epitaxial graphene/Ir(lll) 
system. We show directly that in STM imaging the electronic contribution is pre- 
vailing compared to the topographic one. In the force microscopy and spectroscopy 
experiments we observe a variation of the interaction strength between the tip and 
high-symmetry places within the graphene moire supercell, which determine the ad- 
sorption cites for molecules or metal clusters on graphene/Ir(lll). 
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Graphene layers on metal surfaces have been attracting the attention of scientists since 
several decades, starting from middle of the 60s, when the catalytic properties of the close- 
packed surfaces of transition metals were in the focus of the surface science research [IHl]- 
The demonstration of the fascinating electronic properties of the free-standing graphene [5], 
E], renewed the interest in the graphene /metal systems, which are considered as the main and 
the most perspective way for the large-scale preparation of high-quality graphene layers with 
controllable properties [THTO]. For this purpose single-crystalline as well as polycrystalline 
substrates of 3d — 5d metals can be used. 

One of the particularly exciting questions concerning the graphene/metal interface is the 
origin of the bonding mechanism in such systems [SHU [H]. This graphene-metal puzzle is 
valid for both cases: graphene adsorption on metallic surfaces as well as for the opposite sit- 
uation of the metal deposition on the free-standing or substrate-supported graphene. In the 
latter case the close-packed surfaces of Ad and 5d metals are often used as substrates fl2[ [T3] . 
A graphene layer prepared on such surfaces, i.e. Ru(OOOl) P^l - fTB] . Rh(lll) |T31 [IZl [TH] . 
Ir(lll) [ISl 120], or Pt(lll) [211 122], forms so-called moire structures due to the relatively 
large lattice mismatch between graphene and metal substrates. As a consequence of the 
lattice mismatch the interaction strength between graphene and the metallic substrate 
is spatially modulated leading to the spatially periodic electronic structure. Such lateral 
graphene superlattices are known to exhibit selective absorption for organic molecules [23| 
or metal clusters [21] . Especially, the adsorption of different metals - Ir, Ru, Au, or Pt - on 
graphene/Ir(lll) has been intensively studied showing a preferential nucleation around the 
so-called FCC or HCP high-symmetry positions within the moire unit cell [121 125]. In the 
subsequent works [23 [26] this site-selective adsorption was explained via local sp^ to sp^ 
rehybridization of carbon atoms with the bond formation between graphene and the cluster. 
However, a fully consistent description of the local electronic structure of graphene/Ir(lll), 
the observed imaging contrast in scanning probe experiments and the bonding mechanism 
of molecules or clusters on it is still lacking, motivating the present research. 

Here we present the systematic studies of the graphene /Ir(l 11) system by means of density 
functional theory (DFT) calculations and scanning tunnelling and atomic force microscopy 
(STM and AFM) performed in constant current / constant frequency shift (CC / CFS) 
and constant height (CH) modes. The obtained results for the graphene/Ir(lll) system 
allow to separate the topographic and electronic contributions in the imaging contrast in 
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STM and AFM as well as to shed light on the spatially modulated interaction between 
graphene/Ir(lll) and the metallic STM/AFM tip, which is of paramount importance for 
the understanding of absorption of metals on top of graphene/Ir(lll) as well as of similar 
graphene-metal systems. 

Results 

The unit cell of graphene on Ir(lll) is shown in Fig. [l|a) with the corresponding high- 
symmetry local arrangements of carbon atoms above Ir layers marked in the figure: ATOP 
(circles), FCC (squares), HCP (down-triangles), and BRIDGE (stars). The DFT-D2 
optimized local distances between graphene and Ir(lll) are 3.27A (HCP), 3.28A (FCC), 
3.315 A (BRIDCE), 3.58 A (ATOP). This result is very close to the recently published 
equilibrium structure for this system [20] • The similar distances for HCP, FCC, and 
BRIDGE positions can be related to the fact that in all these cases one of the carbon 
atoms in the graphene unit cell is placed above Ir(S) atom defining the local interaction 
strength for the particular high-symmetry position. The obtained distances are very close to 
those between carbon layers in pure graphite and this was explained by a binding interaction 
dominated by van der Waals effects between graphene and Ir(lll) that is modulated by weak 
bonding interactions at the FCC and HCP places and anti-bonding chemical interaction 
around ATOP positions . As a result a small charge transfer from the graphene tt states 
on Ir empty valence band states is detected in calculations and the Dirac point is shifted by 
^ 100 meV above the Fermi level [Ep) that is close to the data obtained by photoelectron 
spectroscopy [31] [see Fig. ^c) and discussion below] . 

The graphene/Ir(lll) system is a nice example of the moire structure, which is easily 
recognisable in LEED and on the large scale STM images shown in Fig. [l](b). The extracted 
lattice parameter of this structure from LEED and STM is 25.5 A and 25.2 A, respectively, 
that is in good agreement with previously published data [33j- In most cases, for the typical 
bias voltages used in STM imaging, the graphene/Ir(lll) structure is imaged in the so-called 
inverted contrast [Fig. [l|^c,d)] [3S|, when topographically highest ATOP places are imaged 
as dark and topographically lowest FCC and HCP as bright regions. This assignment 
was initially done in Ref. [33], where the STM topography of neighbouring regions of the 
clean and graphene-covered Ir(lll) surface were imaged. In our studies we perform the 



3 



"on-the-fly" switching between CC STM and CFS AFM imaging during scanning [Fig. [T](d)] 
where we observe the inversion of the topographic contrast z{x, y) (see also Fig. SI of the 
supplementary material [55]). This becomes clearly evident around areas marked by arrows 
in Fig. [l|d) where the darkest contrast in CC STM becomes the brightest one in CFS AFM 
for the ATOP position. In the latter case the darkest areas correspond to HCP sites that 
correlates with the calculated height changes between lowest and highest carbon positions. 
The pronounced difference between the FCC and HCP areas can be explained by fact that 
additional bias voltage was applied in order to increase the atomically-resolved contrast. 

The inversion of the imaging contrast was also detected in CC STM images when the 
bias voltage is changed from —0.5 V to —1.8V during scanning [Fig. [2|a)]. The images 
taken at positive voltages are similar to the lower part of the figure Fig. [2]^a) (see Fig. S2 of 
the supplementary material [35]). The simulations of STM images using the experimental 
conditions show very good agreement with the obtained data [Fig. |2](a,b) and Fig. S2]. The 
presented STM results demonstrate the big difference in the local electronic structure for 
high-symmetry positions of carbon atoms on Ir(lll). That can points out on the difference 
in the local adsorption strength for these places. 

In order to prove this assertion we have performed the force microscopy and spectroscopy 
experiments on graphene/Ir(lll). Fig. [s] shows (a,b) the frequency shift of an oscillating 
scanning sensor as a function of a distance from the sample, Af{d), and (c) the corresponding 
tunnelling current, I{d), measured in the unit cell of graphene/Ir(lll) along the path marked 
in the CC STM image shown as an inset of (a). During these measurements the tunnelling 
current was used for the stabilisation of the feedback loop, allowing to determine the relative 
2;-position of the A/ and / curves. Following the sequence: Fz{d) = —dE{d)/dd, Af{d) = 
—fo/2ko ■ dFz{d)/dd, where E{d) and Fz{d) are the interaction energy and the vertical force 
between a tip and the sample, respectively, /o and fco are the resonance frequency and the 
spring constant of the sensor, correspondingly, we can separate repulsive, attractive, and 
long-range electrostatic and van der Waals contributions in AFM imaging. The imaging 
contrast at distances more than 5 A from the surface is defined by long-range van der Waals 
interactions which are insensitive to the atomically resolved structure of the scanning tip 
and the sample. The short-range chemical interactions around the minimum of the A/ curve 
are dominated either by repulsive (left-hand side) or attractive (right-hand side) forces and 
give the atomically-resolved site-selective chemical contrast. 
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The presented ^f{d) curves [Fig. [3]^a,b)] clearly show the site-selective interaction in the 
unit cell of the graphene/Ir(lll) system. First of all, the absolute value of the maximal 
frequency shift is higher for the FCC and HCP areas compared to ATOP ones indicating 
the stronger interaction of graphene with the tungsten STM/AFM tip for the former carbon 
positions. We would like to note, that the same trend will be observed for the adsorption of 
metals on graphene/Ir(lll). Secondly, the difference in the positions of Af{d)-cmve minima 
is 0.097 nm [Fig. ^h)] compared to the value of graphene corrugation of 0.067 nm obtained 
from the CC STM image [inset of Fig. |3](a)], also pointing on the difference in the interaction 
strength for two carbon positions, ATOP and FCC. 

The measured Af{d) curves were used to track and explain the contrast obtained in CH 
AFM experiments (Fig. |4]and Fig. 4 in Supplementary material [35j). In this experiments 
the z-coordinate of the oscillating sensor was fixed during scanning and A/(x, y) and I{x, y) 
maps were collected. As can be clearly seen from the presented results the imaging contrast 
in CH A/ images is changing as a function of 2;-coordinate and it is fully inverted for two limit 
2;-positions [di = — 0.033 nm and d2 = +0.227 nm, see Fig. [3]^b)] used for the imaging (upper 
row in Fig. |4]). In the region of distances between the tip and the graphene/Ir(lll) sample 
where chemical forces determine the interaction, on the right-hand side of the A/ curve the 
interaction is attractive giving the more negative frequency shift for the ATOP positions as 
located closer to the oscillating tip compared to the FCC and HCP areas. On the other 
side, on the left-hand side of the curve the repulsive interaction starts to contribute in the 
interaction giving the smaller frequency shift for ATOP compared to other positions that 
can again be explained by the different distances between tip and the sample for different 
high-symmetry positions as well as by the different interaction strength at these places. At 
the same time the contrast for the current map I{x,y) is always the same, inverted contrast 
(lower row in Fig. |4]). Taking into account that during CH imaging at +50meV in Fig. |4] 
the imaging contrast for I{x,y) is inverted due to the higher DOS for the FCC and HCP 
positions and that variation of the distance for CH imaging and CC STM imaging (when 
Ut is changed from —0.5 V to —1.8 V) is nearly the same, one can separate the topographic 
and electronic contributions into imaging at different biases and distances. 

Here we would also like to note that the change of the bias voltage during CH AFM 
imaging does not lead to any changes in the imaging contrast for A/ (Fig. |5|. The first 
row shows the two small-scanning range atomically resolved Af{x,y) images acquired on 
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the same place of the graphene/Ir(lll) sample with opposite signs for the bias voltage: 
the imaging contrast is the same with the slight variation of the imaging scale that can be 
explained by the small drift of the oscillating tip. The first look on the I{x, y) map (lower 
row) might give an impression that the contrast is fully inverted. However, the absolute 
value of the tunnelling current is the same and only more negative values of the tunnelling 
current are shown as darker areas in the image for current. 

Discussion 

The graphene/Ir(lll) system was studied with DFT methods in oder to deeply under- 
stand the effects of the electronic structure on the adsorption properties of this system. Here 
we present the comparison between experimental and theoretical STM data. The STM im- 
ages are calculated using the Tersoff-Hamann formalism ^B], in its most basic formulation, 
approximating the STM tip by an infinitely small point source [iHl [371 EH] • In these simula- 
tions the constant current condition was fulfilled that leads to the increasing of the distance 
between the sample and a tip from 2.50 A for Ut = —0.5 V to 3.21 A for Ut = —1.8 V 
[Fig. [2|^b)]. Taking into account that at negative bias voltages the electrons tunnel from 
occupied states of the sample into unoccupied states of the tip and that tunnelling current 
is proportional to the surface local density of states (LDOS) at the position of the tip we 
can clearly identify the states in the valence band of graphene/Ir(lll) which are respon- 
sible for the formation of the imaging contrast (see Fig. S3 of the supplementary material 
for the identification of the local valence band states states). [The surface Brillouin zone 
(BZ) of the graphene/Ir(lll) system is approximately 10 times smaller compared to BZ of 
free-standing graphene giving the maximal wave- vector of electrons of A;|| ~ 0.17A~^ for the 
ii'-point, i. e. in this case one has to consider a small region in the reciprocal space around 
the F-point. In this case we can assign the features in the tunnelling current as originating 
from the corresponding peaks in DOS of the whole system (in general, the tunnelling current 
depends on the k\\ of electronic states participating in the tunnelling process and the elec- 
tronic states with a large parallel component of the wave vector give a small contribution to 
the current)]. Therefore, the so-called inverted contrast for the small bias voltages, negative 
or positive, can be explained by two features in LDOS for the carbon atoms in the FCC 
position located at E — Ep = — 0.53eV and E — Ep = +0.32 eV, respectively [Fig. |2|^c)]. 
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The next two maxima in LDOS for C ATOP positions located at E — Ep = — 0.95eV and 
E — Ep = — 1.49eV [Fig. ^c)] explain the fact that these places become brighter for higher 
bias voltages in experimental and simulated CC STM images [Fig. [2]^a,b)]. 

The calculated LDOS for graphene/Ir(lll) compared with DOS for the free-standing 
graphene [(10 x 10) unit cell] is shown in Fig. ^c) together with the electron density 
difference distribution presented as an inset. The obtained DOS results demonstrate the 
slight difference in the position of main features (vr-states at — Ep ^ 6.1 eV, a-states at 
E — Ep ^ 3.2 eV, and the M-point-derived Kohn anomaly at — Ep ^ 2.2 eV) as well as 
the Dirac cone between two systems, that is assigned to the p-doping of graphene on Ir(lll). 
The DOS features responsible for the formation of the STM contrast can be considered as 
a result of the hybridisation between graphene vr states and the valence band states of Ir 
around the particular high-symmetry adsorption positions (see Fig. S3 of the supplementary 
material for the identification of the local valence band states states). Here we would like 
to note, that one can also consider the increased LDOS around Ep for FCC, HCP, and 
BRIDGE positions as a hint indicating that these places can be a nucleation centers for 
adsorbed atoms. 

In oder to get better insight in the results obtained during AFM experiments we per- 
formed simulations of these data in the framework of DFT formalism (we would like to note 
that a qualitative description of the buckled graphene systems was performed in Ref. |39] . 
where the model Lennard- Jones potential was used to model tip-sample interaction). Due to 
the large size of the unit cell of the graphene/Ir(lll) system we selected two approaches. In 
the first one, which requires less computational resources, the tip-sample force is expressed 
as a function of the potential Vts{r) on the tip due to the sample: Fts oc — V[|VVjs(r)p] 
However, this method does not take into account the geometrical and electronic structure of 
the scanning tip and thereby does not allow to get absolute values for the force or frequency 
shift and the correct distance between the tip and the sample, and only qualitative result 
for the attractive region of the interaction can be obtained. The result of such simulations 
for graphene/Ir(lll) is shown in Fig. [6]^a), which is in rather good agreement with exper- 
imental data: the regions around ATOP positions are imaged as dark compared to other 
high-symmetry cites imaged as brighter contrast. Unfortunately, the information about re- 
pulsive region of the forces can not be obtained from such calculations. However, if the 
model corrugation-dependent repulsive potential is added, then the resulting simulated A/ 
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curve reproduces qualitatively the experimentally obtained results. 

For the purpose to reproduce our data in a more quantitative way, the interaction between 
the tip and the graphene/Ir(lll) system was simulated within the second approach, where 
W-tip is approximated by the 5-atom pyramid as shown in Fig. ^h). The results of these 
calculations are compiled in Fig. [6]^c-e). The interaction energy (system was rigid without 
allowing to relax) between model W-tip and the surface was calculated for two limiting 
places of graphene/Ir(lll), ATOP and FCC. The calculated points are shown by filled 
rectangles and circles in Fig. ^c) for the FCC and ATOP positions, respectively. The 
Morse potential was used to fit the calculated data as the most suitable for the graphene- 
metal systems [HI |12]. The resulting curves are shown by solid lines in the same figure. The 
force and the frequency shift curves were calculated on the basis of the obtained interaction 
energy curves according to formulas presented earlier. The results are shown in Fig. [6](d) 
and (e), respectively. 

The interaction energy curves [Fig. ^c)] are very similar for the two high-symmetry ar- 
eas and shifted with respect to each other by 0.40 A reflecting the height difference between 
ATOP and FCC regions. The maximal absolute value of energy is slightly higher for FCC 
by 0.015 eV. However, our calculations do not take into account the relaxation of the system 
during the tip-sample interaction. Therefore the actual values can be different. Calcula- 
tion of the force and the frequency shift via derivation of the respective curves leads to the 
clear discrimination in the value of the extrema as well as to the further separation of their 
positions and for the A/ curves the extrema are separated by 0.45 A. The resulting values 
for the frequency shift for our model W-tip are approximately 2 times larger compared to 
experimental values. But, the shape and the trend for the inversion of the imaging contrast 
are clearly reproducible in our calculated data. The following reasons can explain the dif- 
ference between experiment and theory: (i) small W-cluster modelling the actual W-tip, (ii) 
system was not relaxed during calculations, (iii) the actual atomic structure of the tip during 
experiment is unknown, (iv) small uncertainty in the tip position during measurements of 
A/ curves. However, in spite of these simplifications, the presented calculation results allow 
us to confirm and better understand the data obtained during the AFM measurements. 

In conclusion, we performed the systematic investigations of the geometry, electronic 
structure and their effects on the observed imaging contrast during STM and AFM experi- 
ments on graphene/Ir(lll). Our results obtained via combination of DFT calculations and 



8 



scanning probe microscopy imaging allow to discriminate the topographic and electronic 
contributions in these measurements and to explain the observed contrast features. We 
found that in STM imaging the electronic contribution is prevailing compared to the to- 
pographic one and the inversion of the contrast can be assigned to the particular features 
in the electronic structure of graphene on Ir(lll). Contrast changes observed in constant 
height AFM measurements are analyzed on the basis of the energy, force, and frequency 
shift curves reflecting the interaction of the W-tip with the surface and are attributed to the 
difference in the height and the different interaction strength for high-symmetry cites within 
the moire unit cell of graphene on Ir(lll). The presented findings are of general importance 
for the understanding of the properties of the lattice- mismatched graphene/metal systems 
especially with regard to possible applications as a template for molecules or clusters. 

Methods 

DFT calculations. The crystallographic model of graphene/ Ir (111) presented in Fig. [l|^a) 
was used in DFT calculations, which were carried out using the projector augmented plane 
wave method |27], a plane wave basis set with a maximum kinetic energy of 400 eV and the 
PBE exchange-correlation potential [2H], as implemented in the VASP program [22! • The 
long-range van der Waals interactions were accounted for by means of a semiempirical DFT- 
D2 approach proposed by Grimme |30]. The studied system is modelled using supercell, 
which has a (9 x 9) lateral periodicity and contains one layer of (10 x 10) graphene on 
a four-layer slab of metal atoms. Metallic slab replicas are separated by ca. 20 A in the 
surface normal direction. To avoid interactions between periodic images of the slab, a dipole 
correction is applied |3T]. The surface Brillouin zone is sampled with a (3 x 3 x 1) /c-point 
mesh centered the F point. 

STM and AFM experiments. The STM and AFM measurements were performed in two 
different modes: constant current (CC) or constant frequency shift (CFS) and constant 
height (CH). In the first case the topography of sample, z{x,y), is studied with the corre- 
sponding signal, tunnelling current {It) or frequency shift (A/), used as an input for the 
feedback loop. In the later case the 2;-coordinate of the scanning tip is fixed with a feed- 
back loop switched off that leads to the variation of the distance d between tip and sample. 
In such experiments lT{x,y) and Af{x,y) are measured for different 2;-coordinates. The 
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STM/AFM images were collected with Aarhus SPM 150 equipped with KolibriSensor"'"'^ 
from SPECS [ISl E2] with Nanonis Control system. In all measurements the sharp W- 
tip was used which was cleaned in situ via Ar+-sputtering. In presented STM images the 
tunnelling bias voltage, Ut, is referenced to the sample and the tunnelling current, It, is 
collected by the tip, which is virtually grounded. During the AFM measurements the sensor 
was oscillating with the resonance frequency of /o = 999161Hz and the quality factor of 
Q = 45249. The oscillation amplitude was set to A = 100 pm or A = 300 pm. 

Preparation of graphene/Ir(l 1 1) . The graphene/Ir(lll) system was prepared in ultra- 
high vacuum station for STM/AFM studies according to the recipe described in details 
in Ref. p3] via cracking of ethylene: T = 1100° C, p = 5 x lO^^mbar, t = 5min. This 
procedure leads to the single-domain graphene layer on Ir(lll) of very high quality that was 
verified by means of low-energy electron diffraction (LEED) and STM. The base vacuum 
was better than 8 x 10^^^ mbar during all experiments. All measurements were performed 
at room temperature. 
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Fig. SI. STM/AFM images of the graphene/lr(lll) system obtained via switching 
between CC STM and CFS AFM modes "on-the-fly" . Scanning parameters: (a) STM: 
Ut = +461 meV, It = 7.6 nA, AFM: A/ = -675 mHz (attractive regime), (b) STM: 
Ut = -201 meV, It = 7.6 nA, AFM: A/ = +675 mHz (repulsive regime). 
Fig. S2. (a) STM image of graphene/Ir(lll) obtained at Ut = +0.7 V and It = 15 nA. 
(b) Calculated STM image obtained in the framework of the Tersoff-Hamann formalism via 
integration of the valence band states in the range E — Ep = 0... + 0.5eV. The corresponding 
distance between tip and the sample surface is 2.2 A. 

Fig. S3. Band structures of graphene/Ir(lll) calculated for the corresponding arrange- 
ments of the expanded graphene layer on Ir(lll) in (1 x 1) structure. The distance 
between graphene and Ir(lll) is fixed to the corresponding spacing in the graphene/Ir(lll) 
nanomesh. The thickness of the energy bands corresponds to the higher atom- and 

state-projected contribution to the corresponding band. 

Fig. S4. The sequence of the CH AFM images (A/ and /) obtained at different 2;-positions 
of the scanning tip and correspondingly at different distances d between tip and the surface. 
The bias voltage apphed to the oscillating tip during scanning is +50 meV. 
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FIG. 1: (a) Crystallographic structure of (10 x 10) graphene on (9 x 9) Ir(lll). The high- 
symmetry places are marked by circle, rectangle, down-triangle, and stars for ATOP, FCC, HCP, 
and BRIDGE positions, (b) Large scale (162 x 55 nm^) STM image of graphene/Ir(lll). The 
inset shows the corresponding LEED image obtained at 71 eV. (c) Atomically resolved STM image 
of graphene/Ir(lll) showing the inverted contrast, (d) A combined STM/AFM image with a 
switching between CC STM and CFS AFM imaging during scanning. The scan sizes are 5.2 x 
5.2 nm2 {Ut = +550 meV, It = 540 pA) and 10 x 10 nm^ (STM: Ut = +461 meV, It = 7.6 nA, 
AFM: A/ = — 675mHz), respectively. 
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FIG. 2: Experimental (a) and calculated (b) CC STM images of graphene/Ir(lll) obtained at 
—0.5 V (bottom) and —1.8 V (top) of the bias voltage (It = 15 nA). (c) Carbon-site projected 
density of states calculated for all high-symmetry places of graphene/Ir(lll) (zoom around Ep 
is shown as an upper inset). Difference electron density, Ap{r) = Pgr/ir{iii){''') ~ Piri^) — Pgr{r), 
plotted in units of e/A^ is shown as an inset. 
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FIG. 3: Frequency shift (a,b) of the oscillating KolibriSensor and the corresponding tunnelling 
current (c) as a function of the relative distance d between tip and the graphene/Ir(lll) sample. 
These curves were acquired with respect to the set-point for the CC STM imaging. Inset of (a) 
shows the corresponding STM image (Ut = +50 mV, It = 400 pA) with the path where A/ and 
/ data were measured. 
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FIG. 4: Constant height images of graphene/Ir(lll), A/(x, y) (upper row) and /(x, y) (lower row), 
obtained at two different heights di and d2 of the KohbriSensor above the surface (see Fig. [3]). The 
size of the scanning area is 10.5 x 10.5 nni'^. 
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t/r= + 400 mV Ut= -400 mV 

-1.00Hz -1.19Hz 




277 pA -385 pA 



FIG. 5: Atomically-resolved constant height images of graphene/Ir(lll), Af{x,y) (upper row) 
and I{x, y) (lower row), coUected at opposite signs bias voltages. The size of the scanning area is 
4.3 X 4.3 nm^. 
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FIG. 6: (a) Simulated CH AFM image of the graphene/Ir(lll) system according to the ap- 
proach suggested in Ref. llOj. (b) Schematic representation of the geometry of the 5-atom-W- 
tip/graphene/Ir(lll) system used in simulation of the frequency shift curves. Interaction energy 
(c) and the force (d) between 5-atom W-tip and graphene/Ir(lll) calculated for two limit places 
ATOP and FCC. (e) Frequency shift as a function of distance between the model 5-atom W-tip 
and graphene/Ir(lll). 
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Fig. SI. STM/AFM images of tlie graphene/Ir(lll) system obtained via switching be- 
tween CC STM and CPS AFM modes "on-t he-fly" . Scanning parameters: (a) STM: 
Ut = +461 meV, It = 7.6 nA, AFM: A/ = -675 mHz (attractive regime), (b) STM: 
Ut = — 201meV, It = 7.6 nA, AFM: A/ = +675 mHz (repulsive regime). 
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Fig. S2. (a) STM image of grapliene/Ir(lll) obtained at Ut = +0.7 V and It = 15 nA. 
(b) Calculated STM image obtained in the framework of the Tersoff-Hamann formalism via 
integration of the valence band states in the range E — Ep = 0... + 0.5eV. The corresponding 
distance between tip and the sample surface is 2.2 A. 




Fig. S3. Band structures of graphene/Ir(lll) calculated for the corresponding arrange- 
ments of the expanded graphene layer on Ir(lll) in (1 x 1) structure. The distance be- 
tween graphene and Ir(lll) is fixed to the corresponding spacing in the graphene/lr(lll) 
nanomesh. The thickness of the energy bands corresponds to the higher atom- and state- 
projected contribution to the corresponding band. 
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Fig. S4. The sequence of the CH ARVE images [A/ (left column) and / (right column)] 
obtained at different z-positions of the scanning tip and correspondingly at different distances 
d between tip and the surface. The bias voltage applied to the oscillating tip during scanning 
is +50 meV. 
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Fig. S4. (continued) The sequence of the CH AFM images [A/ (left column) and / (right 
column)] obtained at different z-positions of the scanning tip and correspondingly at different 
distances d between tip and the surface. The bias voltage applied to the oscillating tip during 
scanning is +50 meV. 
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